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Abstract 

This  paper  explores  low  observability  flight  path  planning  of  unmanned  air  vehicles  in  the 
presence  of  radar  detection  systems.  The  probability  of  detection  model  of  an  aircraft  near 
an  enemy  radar  depends  on  aircraft  attitude,  range,  and  configuration.  A  detection  model 
is  coupled  with  a  simplified  aircraft  dynamics  model.  The  Nonlinear  Trajectory  Generation 
(NTG)  software  package  developed  at  Caltech  is  used.  NTG  algorithm  is  a  gradient  descent 
optimization  method  that  combines  three  technologies:  B-splines,  output  space  collocation  and 
nonlinear  optimization  tools.  Implementations  are  formulated  with  temporeJ  constraints  that 
eJIow  i>eriods  of  high  observabihty  interspersed  with  periods  of  low  observability.  Illustrative 
examples  of  optimized  routes  for  low  observability  are  presented. 

Keywords:  Trajectory  Generation,  Path  Planning,  Low-Observable,  UAV,  RADAR,  Probability  of  Detection 


1  Introduction 

As  the  development  of  new  unmaimed  vehicles  progresses,  research  focus  on  coordinating  large 
teams  of  these  vehicles  has  highlighted  many  new  challenges  for  control.  Due  to  complexity,  the 
approach  taken  in  the  DARPA  MICA  program  is  to  decompose  the  control  design  in  a  hierarchical 
manner.  The  cooperative  path  planning  layer  routes  the  vehicles.  At  this  level  of  the  controller,  the 
design  goals  are  to  command  the  vehicles  to  pursue  dynamically  feasible  routes,  such  as  speed  and 
turning  rate  limits,  while  avoiding  threats  and  collisions  with  other  vehicles.  Driven  by  battlespace 
management  needs,  real-world  threats  can  be  much  more  complex  than  those  modeled  in  prior 
trajectory  generation  work.  The  detectability  of  an  aircraft  traveling  neeir  an  enemy  radar  depends 
on  more  than  just  the  distance  to  the  radar;  it  depends  on  the  aircraft  attitude  and  configuration 
as  well.  This  feature  in  the  threat  model  introduces  path  dependencies  as  well  as  sharp  gradients 
into  the  underlying  optimization  problems,  and  it  presents  new  challenges  for  trajectory  generation 
techniques.  In  this  paper,  we  investigate  the  use  of  the  nonlineeir  trajectory  generation  (NTG) 
method  as  a  solution  to  the  low-observability  path  planning  challenge. 

For  the  low-observability  problem,  the  use  of  NTG  is  motivated  by  recent  extensions  in  its 
ability  to  deal  with  temporal  constraints.  In  [5],  NTG  is  used  for  a  missile  intercept  problem.  In 
[4],  NTG  is  extended  to  a  multi-vehicle  problem  with  precedence  constraints  such  as  look  after  strike 
and  simultaneous  strike.  The  low-observability  routing  problem  could  be  considered  to  be  temporal 
in  nature,  by  allowing  periods  of  high  observability  interspersed  with  periods  of  low  observability. 
This  is  desirable  because  of  the  way  the  enemy  fire-control  systems  work.  Although  it  might  not 
be  possible  to  get  close  to  the  enemy  targets  while  maintaining  low-observability  at  aU  times;  by 
strategically  flying  low-observable  paths  for  part  of  the  time,  it  may  be  possible  to  drive  the  enemy 
systems  into  a  condition  called  lock-loss.  This  condition  aborts  the  enemy  fire  control  plans  after  a 
specified  time  of  no  detection. 

Computational  efficiency  as  well  as  the  capability  to  enforce  more  reeilistic  constraints,  are  two 
additional  motivations  that  have  prompted  the  choice  for  NTG  in  past  problems.  The  method 
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combines  three  technologies-  differential  flatness,  splines  and  nonlinear  optimization  and  has  been 
extensively  investigated  [6,  7,  3]  .  Differentially  flat  systems  have  the  property  that  the  input 
vEiriables  and  states  can  be  written  in  terms  of  the  output  variables  and  their  derivatives.  This 
can  aid  in  computationally  efficiency  by  eliminating  the  need  to  integrate  as  in  shooting  methods 
[1].  The  use  of  splines  also  addresses  computational  issues  by  eJlowing  complicated  functions  to  be 
written  with  low  order  polynomials  that  are  active  over  distinctive  intervals.The  NPSOL  [8]  package, 
which  uses  sequential  quadratic  programming,  is  used  in  this  work  as  a  nonlinear  optimization  tool. 

In  this  paper,  we  investigate  the  use  of  NTG  for  low-observability  trajectory  generation  for 
unmanned  air  vehicles.  The  aircraft  model  and  detection  models  are  described  in  the  second 
section.  The  third  section  presents  an  NTG  approach  formulated  for  low  observability.  The  fourth 
section  contains  examples.  The  exeunples  demonstrate  the  ability  of  NTG  to  converge  to  solutions 
which  constredn  observability  to  acceptable  levels;  however,  they  also  highlight  some  difficulties 
with  the  approach. 

2  Model 

This  section  presents  the  model.  These  models  are  based  on  those  developed  as  part  of  the  DARPA 
MICA  progreun  by  [2].  In  this  progrsun,  as  the  problems  become  progressively  more  challenging, 
the  models  become  more  refined  to  capture  more  realistic  features.  For  the  purposes  of  this  report, 
the  models  used  here  are  simplified. 

The  two  main  components  are  the  aircraft  and  the  detection  models,  as  shown  in  figure  2.  The 
aircraft  model  shown  is  a  simple  point  mass  model.  Throughout  this  paper,  we  assume  that  the 
aircraft  maintains  a  fixed  altitude  and  the  radar  is  on  the  ground.  The  inputs  to  the  model  are  a 
series  of  waypoints  as  well  as  the  average  speeds  required  to  travel  between  waypoints.  Changes 
in  speed  and  heading  occur  instantaneously  at  the  waypoints.  The  output  of  the  aircraft  model  is 
the  aircraft  position  emd  attitude  in  inertial  coordinates  (East,  North,  Up).  Figure  1  illustrates  an 
aircraft  traveling  between  waypoints.  The  equations  below  represent  the  state  equations  for  the 
vehicle  traveling  fi-om  waypoint  i  to  waypoint  j. 


iac  =  UijSin{il>ij)  (1) 

Vac  =  UijCOs(lpij)  (2) 

(3) 

where  Uij  is  the  speed  emd  tp  is  the  heading,  the  angle  between  the  nose  and  north.  Equivalently, 
assuming  that  the  velocity  vector  is  aligned  with  the  nose. 


xp  =  tan  (4) 

Vac 

The  aircraft  position  {xac,  Vaci  Zac)  is  combined  with  the  radar  position  (xr,  yR,  zr)  to  form  a 
vector  from  the  aircraft  to  the  radar. 

~  {.^R  ~  Xac)Xeast  d"  {VR  ~  yac)ynorth  "b  i^R  ~  ^ac)Zyip  (5) 

This  vector  is  then  transformed  to  body  axes.  Assuming  zero  bank  and  pitch  angle,  for  this 
model,  the  transformation  is 


^Rac 

cos{ip) 

sin{ip) 

0  ■ 

{XR  -  Xac) 

Rhody  — 

VRoc 

= 

—sin{ip) 

cos{ip) 

0 

(VR  -  Vac) 

ZRac 

0 

0 

1 

{zR  -  Zac) 

The  inputs  to  the  detection  model  are  the  azimuth,  elevation  and  slant  range  (az,el,Rg).  To 
obtain  these  values  the  vector  Rbody  is  transformed  to  spherical  coordinates  as  follows 
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Figure  1:  The  sdrcreift  model  is  a  point  mass  model. 

az  =  fon-i  ( 

(7) 

— 

\XRacJ 

— 

'(ji,) 

(8) 

=  v/^Lc  +  !4c  + 

(9) 

where  Rg  =  is  the  ground  range.  Note  that  since  heading  is  a  function  of  velocities' 

as  in  equation  4,  then  these  equations  can  be  rewritten  as  functions  of  the  position  and  velocity 

XyX,y,y 

For  the  detection  model,  two  tabl^s-iBe  used,  these  tables  are  based  on  the  small  uav  and 
Long  SAM  radar  model  parameters  ib  [2\J!ae  first  table  computes  the  radar  signature  given  the 
azimuth  and  elevation  in  degrees.  ThelSigimure  is  a  unitless,  intermediary  variable  that  is  related 
to  radar  cross  section.  Note  for  azimuths  within  +/-  30  degrees  (“nose-in”  flight),  the  signature 
values  are  lower  than  for  azimuths  outside  this  range  (“nose-out”  flight).  The  second  table  relates 
the  probability  of  detection  with  the  signature  and  the  slant  range. 

Table  1.  Signature  Values  Respect  to  Azimuth  and  Elevation 


el  /az 

0.0 

-I-/-30 

-I-/-31 

-I-/-180 

0.0 

1.5 

1.5 

5.5 

5.5 

-t-/-20 

2.5 

2.5 

5.5 

5.5 

-I-/-45 

3.5 

3.5 

6.0 

6.0 

-I-/-90 

6.5 

6.5 

6.5 

6.5 

Table  2.  Probability  of  Detection  (Pd)  Values  Respect  to  Range  smd  Signature 
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Figure  2:  The  aircrsrft  dynamics  and  the  detection  model  comprise  two  main  components  of  the 
model. 


Signature 

Pd=.99 

Pd=.5 

Pd=.l 

Pd=.01 

7 

380.1 

481.2 

555.6 

656.6 

6 

213.7 

270.6 

312.5 

369.2 

5 

120.2 

152.2 

175.7 

207.6 

4 

67.6 

85.6 

98.8 

116.8 

3 

38.0 

48.1 

55.6 

65.7 

2 

21.4 

27.1 

31.2 

36.9 

1 

12.0 

15.2 

17.6 

20.8 

At  fixed  altitude,  it  is  interesting  to  consider  how  the  probability  of  detection  varies  with 
ground  range  as  in  figure  {sj.  The  figure  shows  that  “nose  in”  flight  is  signficantly  less  detectable 
than  “nose  out  flight”.  Near  50  km,  the  approximate  range  of  a  homing  missile,  the  detection 
probabilities  vriry  approximately  from  0.4  to  0.65.  Outside  the  80  km  radar  range,  the  detection 
probabilities  are  zero. 

To  aid  in  understanding  how  the  probability  of  detection  varies  along  various  trajectories,  we 
focus  on  an  exEunple  where  the  radar  is  located  at  the  origin  of  an  inertied  coordinate  freune  with 
coordinates  (x  (east),y  (north),z  (up)).  An  aircraft  starts  at  position  (-100  km,  -100  km,  12  km) 
and  may  travel  along  diflferent  flight  paths.  Two  flight  paths  are  shown  in  figure  4.  Positions 
are  shown  in  the  (x  (east),y  (north))  plane.  Times  are  not  indicated  on  the  plot.  The  trajectory 
marked  with  circles  is  a  strmght  line  approach  to  the  origin.  The  heading  is  initially  45  degrees, 
and  the  azimuth  is  always  near  zero.  As  an  airplane  following  the  “circle”  trajectory  flies  closer  to 
the  origin,  the  elevation  angle  has  a  larger  magnitude  and  the  probability  of  detection  rises.  The 
’plus’  trajectory  consists  of  several  segments  with  portions  in  easterly,  northerly  and  straight  line 
approaches  towards  the  radar.  Note  that  the  occasional  sharp  changes  in  the  “plus”  trajectory  occur 
because  the  azimuth  changes  in  and  out  of  the  +/-  30  degree  bound.  In  the  first  table,  azimuths 
outside  this  range  yielded  higher  signatures  and  consequently  higher  probability  of  detections. 

From  the  two  tables,  an  anedytical  observability  model,  generically  shown  in  equation  (10)  can 
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Normal  A/C  Configuration 


Figure  3:  At  fixed  altitude  the  probability  of  detection  is  smaller  for  “nose  in”  flight  with  azimuths 
within  30  deg  and  decreases  with  increasing  ground  range 


be  developed  any  number  of  function  approximation  techniques. 

sig  =  sig(az,  el) 

Pd=Pd(R,sig) 

We  have  used  the  B-spline  curve  fits  because  of  their  flexibility  and  ease  of  computing  their  deriva^ 
tives  to  find  analytical  models  for  the  signature  and  probability  of  detection  data  given  in  tables  (1) 
and  (2),  respectively.  Figure  (5)  illustrates  the  signature  values  given  in  table  (1)  in  3-d.  Figure 
(6)  and  (7)  show  the  result  of  the  fit  function  model  of  the  signature  values  by  B-spline  tensor 
product  functions  amd  the  error  between  the  fit  function  and  actuad  data  values,  respectively. 
Similarly,  Figure  (8)  illustrates  the  probability  of  detection  values  respect  to  rauige  and  signature 
vaJues  given  in  table  (1)  in  3-d.  Figure  (??)  amd  (??)  show  the  result  of  the  fit  function  model  of 
the  probability  of  detection  values  by  B-spline  tensor  product  functions  amd  the  error  between  the 
fit  function  amd  actual  data  values,  respectively. 

Then  the  problem  is  to  find  trajectories  for  the  dynamic  system  in  equation  3  such  that  the 
trajectories  are  optimized  with  respect  to  the  probability  of  detection  models  described  in  this  sec¬ 
tion.  This  section  illustrates  that  even  for  highly  simplified  models,  the  optimization  of  trajectories 
for  low  observability  flight  cam  be  quite  challenging. 

3  Approach 

In  this  section,  we  first  briefly  outline  the  NTG  adgorithm  with  temporal  constraints  and  then 
describe  our  approach  to  the  low  observability  trajectory  optimization  problem  formulation.  Illus¬ 
trative  exaimples  will  be  given  in  the  next  section. 


3.1  Brief  Summary  of  NTG  Algorithm 

The  baseline  NTG  algorithm  hais  been  described  extensively  in  the  Uterature  [5,  7,  6,  3],  therefore 
in  this  section  we  outline  it  briefly. 

The  NTG  software  package  is  based  on  a  combination  of  nonlinear  control  theory,  spline  theory 
amd  sequential  quadratic  programming.  There  are  three  steps  in^^TG  algorithm.  The  first  step 
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Figure  4:  Probability  of  Detection  exhibits  sharp  gradients  and  path  dependencies. 


is  to  use  the  differential  flatness  property  [?]  to  find  a  new  set  of  outputs  of  a  system  so  that  the 
s}T5tem  dynamics  can  be  mapped  down  to  a  lower  dimensional  space.  All  the  states  and  controls  of 
the  system  must  be  recovered  from  the  new  lower  dimensional  representation  of  the  system.  Let’s 
consider  a  nonlinear  system  given  by: 


x  =  f{x,u)  (10) 

where  i  e  iZ"  and  u  €  fT"  represents  the  states  and  the  inputs,  respectively.  It  is  assumed  that 
all  vector  fields  and  fimctions  are  real-analytic.  The  goal  is  to  find  a  minimizing  trajectory  of  the 
system  given  in  Equation  (10)  for  the  following  cost  function  ( J)which  consists  of  initial  (^),  final 
{<f>f)  and  trajectory  (L)  cost  functions; 

J{x,u)  =  <t>o{x(to),u{tfi))-V(i>f{x{tf),u{tf))+  f  '  L{x(t),u(t))dt  (11) 

Jto 

The  s)^tem  given  in  Equation  (10)  is  subject  to  the  following  constraints; 

lbQ<i’o(xito),u{t(i))<ubo  JVb  initial  constraints 

Ibf  <  il)f{x(tf),u{tf))  <  xibf  Nf  final  constraints  (12) 

^  ^  trajectory  constreiints 


The  system  in  Equation  (10)  can  be  mapped  to  a  lower  dimensional  space,  in  which  it  will  be 
easier  £ind  computationally  more  efficient  to  solve  the  optimization  problem,  by  finding  an  output 
z  =  Zi,  -  ■  •  ,Zq  of  the  form 


z  =  A(x,u,u^^\---,u^’'^^  (13) 

where  denotes  ith  derivative  of  u  respect  to  time. 

If  Elquation  (10)  is  differentially  flat  then  the  states  and  inputs  of  the  system,  {x,u),  can  be 
completely  established  from  Equation  (14).  If  no  flat  output  exists  or  one  ceuinot  find  a  flat 

output,  then  (x,  u)  can  still  be  completely  determined  from  the  lowest  dimensional  space  possible 
given  in  Equation  (15).  A  necessary  condition  for  the  existance  of  such  outputs  is  given  in  [?]. 
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(14) 


(x,u)  =  B 

(x,u)  =  (15) 

(x,  u)  =  B2  (z,  •  •  • , 


where  z^*^  denotes  ith  derivative  of  z  respect  to  time. 

The  second  step  in  NTG  is  to  further  represent  these  outputs  in  terms  of  the  B-spline  functions  as 


Pi 

(t)Cf  for  the  knot  sequence  tj,  j  =  ,q 

i=l 


where  Bi^kj{t)  represents  the  B-spline  basis  function  for  the  output  zj  with  the  degree  of  sphne 
polynomial  kj.  €{  represents  the  coefficients  of  the  B-splines,  Ij  is  the  number  of  knot  intervals, 
and  rrij  is  the  number  of  smoothness  conditions  at  the  knot  points,  pj  is  the  number  of  coefficients 
of  the  each  output  given  by 

Pj  =  lj(kj  -  mj)  +  rrij 

Finally  to  solve  the  coefficients  of  the  B-spline  functions  by  sequentied  quadratic  programming 
package  NPSOL,  the  cost  function  and  constraints  given  in  equations  (11)  and  (12),  respectively, 
are  re-formulated  in  terms  of  the  B-spline  coefficients.  Therefore,  the  problem  now  can  be  stated 
as  the  following  nonlinear  programming  form: 

minygijw  F(y)  subject  to  lb  <  c(y)  <  vb 


where  J(x,u) F(y)  and  {<^(x(to),u(to)),  4>/{x{tf),u{tf)),  (^t(x(t),«(t))} -»  c(y) 
y  =  and  M  =  TLiPi 

3.2  Using  Temporal  Constraints  with  NTG 

While  the  NTG  formulation  allows  smy  spatial  constrmnt  to  be  easily  coded  into  the  constraint  set, 
including  tempored  constraints  requires  more  car^[5,  4].  A  key  idea  in  the  temporal  formulation 
of  NTG  is  to  let  event  times  become  state  variable  in  the  optimization.  For  example  suppose  that 
n  events  occur  at  Tf  for  i  =  0, 1,  •  •  •  ,n.  This  allows  precedence  constraints  to  be  included  in  the 
optimization  by  equations  relating  the  event  times.  In  equation  16,  we  define  r  as  a  scaled  time 
vmiable.  It  is  equal  zero  when  time,  t,  is  equal  to  zero,  but  it  is  equed  to  one  when  all  events 
have  occurred  at  T  =  In  the  setup  of  the  optimization  problem,  which  will  be  detailed 

below,  scaled  time  r  rather  than  time  t  is  used.  In  the  optimization  details  below,  the  use  of  new 
time  vmiable,  r,  has  implications  for  the  way  integrals  and  derivatives  are  written.  For  example, 
derivatives  with  respect  to  t  become  derivatives  with  respect  to  r  with  the  chain  rule  as  in  17. 


(16) 
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d  d  dr  _  \  d 

dt  dr  dt  T  dr 


(17) 


As  a  result,  after  introducing  new  state  variable  T,  cost  and  constraint  functions  given  in  equations 
(11)  and  (12)  become 

J(i,u,T)  =  <^(x(0),u(0),T) +  (;6/(i(l),u(l),T)  +  /  L(x(r),u(T),T)dT  (18) 

JO 

lf>0  <  V’o(a:(0),u(0),T)  <  uho  Nq  initial  constraints 

Ibf  <  ipf{x(l),u{l),T)  <  ubf  Nf  final  constraints  (19) 

Ibt  <  rl>t{x{T),u{T),T)  <  ubt  Nt  trajectory  constraints 

There  will  be  also  additional  temporal  constraints  which  can  be  expressed  as  a  set  of  inequalities 
given  below 

Ibr  <  i/>r(T)  <  ubr  Nt  temporal  constraints 


3.3  Formulation  of  the  Low  Observability  Problem 

Now  that  we  have  outlined  the  general  methods,  we  focus  on  how  to  apply  these  methods  to  the 
low  observability  problem.  Define  a  set  of  event  times  T2i  and  T2i+i  for  %  =  0, 1,  •  •  • ,  n  and  the 
number  of  events  2n  +  1  is  specified.  We  will  set  up  events  such  that  periods  of  low  observability 
T2i-i  <t<T2i  are  interspersed  with  times  of  high  observability,  Ta,-  <t<  T2i+i. 

The  system  dynamics  for  this  problem  consist  of  the  vehicle  dynamics  in  equation  (3)  together 
with  the  following  d}rnamics  on  the  new  state  variables 


dT2i+i 

dr 


dT2i 

dr 


=  0 


=  0 


i  =  0,  ...n 


Note  that  the  system  is  differentially  flat.  All  variables  of  interest  can  be  written  as  a  function 
of  the  output  variables,  x,y,T2uT2i+i  and  their  derivatives.  The  observability  model  in  equation 
(10)  is  a  function  of  the  flat  variables,  and  is  used  in  the  optimization  below. 

Next,  we  develop  a  set  of  cost  and  constraint  functions.  Note  that  we  will  develop  a  general 
set,  with  design  options  to  eliminate  some  components  of  the  cost  and  constraint  functions.  A 
cost  function,  J,  is  shown  in  equation  (20).  The  first  term  in  the  cost  function  is  standard  in 
the  temporal  approach.  Since  T  is  the  mission  time  as  shown  in  equation  (21),  this  term  acts  to 
minimize  the  toted  mission  time.  The  second  term  is  a  speed  peneJty.  It  has  edso  been  used  in  the 
past  as  a  way  to  penalize  control  action.  The  third  and  fourth  terms  are  cumulative  penalities  on 
observability.  Since  probability  of  detection  eind  signature  are  always  non-negative,  these  terms  are 
appropriate  and  may  be  useful.  Wu,  Wp  and  Wg  represent  the  weight  functions  on  the  speed, 
probability  and  signature  penalties,  respectively.  The  integral  in  equation  (20)  is  respect  to  the 
scaled  time,  r,  and  has  bounds  from  zero  to  one.  The  cost  function  presented  here  is  a  classic 
tradeoff  between  various  perfomance  measures  £ind  control  action  measures.  In  the  NTG  code 
setup,  the  cost  function  here  would  use  two  subroutines:  wc/,  the  “unintegrated  c(»t  function” 
subroutine  is  used  for  the  costs  that  are  integrated  eind  icf,  the  “initied  cost  function”  may  be  used 
to  add  constant  terms  to  the  cost. 
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where 


T  =  ‘^(T2.  +  T2i+i)  (21) 

t=0 

Next  we  present  a  set  of  constraints,  and  we  divide  the  constraints  into  subsets  based  on  the  type 
of  constradnt  used.  First  there  are  constraints  on  the  initial  conditions  for  the  variables.  There  sue 
hnear  constraints  as  in  equation  (22)  on  the  initial  position  and  event  times.  The  initial  condition 
constraints  on  velocities  become  nonlinear  because  of  the  diflFerentiation  with  respect  to  scaled  time 
shown  above  in  equation  (17).  The  constraints  on  velocities  are  presented  in  equation  (23).  In 
the  NTG  code,  the  linear  initial  conditions  axe  set  up  with  a  variable,  while  the  nonhnear  initial 
conditions  are  written  with  a  subroutine,  nlicf. 


a:(0)  <  x(T)|r=o  <  a;(0) 

1/(0)  <  y(T)  |r=o  <  y(0) 


minT2i  <  T2,(t=o  <  maxT2i 


minT2i+\  <  <  m<ixT2i+i  (22) 

rx(0)  <  —  |r=o  <  Tx(0) 

Ty(0)<— |.=o<ry(0)  (23) 


where  (•)  =  ^ 

Speed  and  radius  of  curvature  limits  are  nonlinear  functions  of  scaled  time.  Constraints  for  these 
are  shown  in  equations  (24)  and  (25).  Note  that  the  equations  are  scaled  to  aid  in  convergence  of 
the  nonhnear  optimization  codes.  These  constradnts  are  written  with  the  subroutine  nltcf  in  the 
NTG  code. 


(24J 


dx  <i^y  __  dy  d?x 

^  3r  dr^  3r  ^  j  f25) 

Finally,  constraints  on  observability  can  be  included.  Consider  a  constreunt  on  signature  as 
in  equations  (26).  In  order  to  circumvent  difQculties  with  discontinuous  constreiints,  which  are 
not  permitted  with  the  NTG  method,  a  single  constraint  on  signature  is  represented  with  two 
continuous  constraints.  During  the  different  low  or  high  observability  events,  only  one  of  these 
two  constraints  is  active.  During  low-observable  times  the  signature  is  constrained  to  be  low  as 
in  the  first  equation  in  the  set  of  equations  in  (26).  The  second  constraint  in  the  second  line  of 
(26)  is  present,  but  is  automatically  satisfied  if  the  first  constraint  is  satisfied.  During  the  times 
of  high-observability,  the  fourth  equation  in  (26)  allows  the  signature  to  reach  its  high  value.  The 
constraint  in  the  third  line  of  (26)  is  present,  but  is  automatically  satisfied. 


T2t  <t  <  T2i+i  0  <  sig  <  sigi 

0  <  sig  <  sign 

T2i+i  <t<  T2i+2  0  <  sig(T2i)  <  sigi 

0  <  siy  <  sign  (26) 
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Constraints  on  probability  of  detection  can  be  included  in  a  similar  fashion  shown  in  equation 
(27). 

32t  <  t  <  T^i+x  0  <  Pd  <  Pd^ 

0<Pd<P<ij, 

T2i+X  <t  <  T2i+2  0  <  Pd(T2i)  <  PdL 

0<Pd<Pd„  (27) 

Note  that  the  constraints  on  signature  and  probability  of  detection  are  nonlinear  function  of 
scaled  time  and  are  represented  in  the  NTG  code  simileirly  to  speed  and  curvature  limits. 

The  last  constraint  in  equation  (28)  ensures  that  the  vehicle  reaches  the  destination  waypoint 
within  a  radius  bound  Cmin  and  emax-  This  constraint  is  nonlinear,  and  it  is  a  function  of  the  final 
scaled  time,  r  =  1,  which  corresponds  to  unknown  actual  final  time  oft/. 

g2 

^  (^(1)  “■  ^destinationwaypointi^))  +  (y(l)  ~  ydestinationwaypointi^))  ^  1  (28) 

^max 

4  Examples 

In  this  secton,  illustrative  examples  to  our  approach  for  the  real  time  low  observabUity  nonlinear 
trajectory  generation  problem  are  given.  We  assume  that  altitude  is  fixed  at  12  km.  The  initial 
location  of  the  aircraft  is  at  (-100  km,  -100  km),  and  there  is  a  radar  located  at  the  origin.  Note  that 
the  difficulties  with  flying  trajectories  in  these  examples  were  highlighted  earlier  in  figure  (4).  In  the 
examples  below,  the  aircraft  travels  to  a  target  located  at  (5  km,  -15  km).  The  examples  highlight 
the  ability  of  the  temporal  methods  to  design  for  distinct  periods  of  high  and  low  observability. 
However,  they  also  demonstrate  sensitivity  to  the  initial  coarse  route. 

Figures  (9), (10),  and  (11)  result  from  a  simulation  in  which  the  initial  coarse  route  input  to  the 
NTG  ailgorithm  is  at  (-100km,  -15  ki|i)^^ure  (9)  shows  that  the  optimized  path  first  circles  through 
the  left  side  of  the  radeir,  thereby  constraining  the  signature  value  of  the  aircraft  during  the  low- 
observability  event  as  shown  in  Figme  (10).  In  this  example  the  low  signature  boimd  is  at  4.2*. 
After  ?sec,  the  signature  is  allowed  to  increase  to  the  maximum  value  at  6.5.  The  azimuth  angle 
of  the  corresponding  flight  path  is  given  in  Figure  (11). 

If  NTG  is  initially  fed  with  a  better  course  route  (0,  -100  km),  then  NTG  prodjmes  a  different 
solution  as  illustrated  in  figures  (12), (13),  and  (14).  In  this  example  the  UAV  flies^ose  in”,  with 
low  azimuths,  directly  toward  the  teirget.  The  signature  values  are  lower  than  in  the  first  example. 

Similar  example  set  also  included  in  Figures  (15)  through  (22).  The  examples  in  figures  (21)  through 
(22)  highlight  that  the  feature  of  the  temporal  method  where  the  signature  is  allowed  to  be  high 
for  a  period  of  time. 

5  Conclusions 

We  will  work  on  this  section 

Because  NTG  is  a  gradient  descent  method,  solutions  that  axe  local  minima  are  possible.  Future 
work  is  investigating  the  use  of  NTG  combined  with  global  approaches  to  address  these  issues. 

Research  supported  as  part  of  the  MICA  program  for  DARPA  and  AFRL,  under  contract 
insertcontractnumber.  ????? 
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actual  aignature  value 


Figure  5;  The  Signature  Data  Versus  Elevation  and  Azimuth  Angles. 
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model  signature  value 


Aulo  H  wMi  tensor  piDduci  apHnes 


Figure  6; 


Eiror  batiMsn  aduel  data  and  ^lin*  fit 


Figure  7;  The  Error  Between  Signature  Data  and  Fit  Model. 
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Figure  8:  The  Probability  of  Detection  Data  Versus  Range  and  Signature  Values. 
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Ground  Track 


Figure  9:  Ground  Track  of  the  Optimized  Path  for  Low  Observability 


Figure  10:  Signature  Value  Across  the  Optimized  Path  for  Low  Observability 
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Azimulh 


Figure  11:  Azimuth  Angle  Across  the  Optimized  Path  for  Low  Observability 


Ground  Track 


Figure  12:  Ground  Track  of  the  Optimized  Path  for  Low  Observability  with  Better  Initial  Course 
Route 
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Si 


Figure  13:  Signature  Value  Across  the  Optimized  Path  for  Low  Observability  with  Better  Initial 
Course  Route 


Figure  14;  Azimuth  Angle  Across  the  Optimized  Path  for  Low  Observability  with  Better  Initial 
Course  Route 


18 


Ground  Track 


Figure  15:  Ground  T^ack  of  the  Optimized  Path  for  Low  Observability 


Figure  16:  Signature  Value  Across  the  Optimized  Path  for  Low  Observability 
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Figure  19:  Signature  Value  Across  the  Optimized  Path  for  Low  Observability 


Azimuth 


Figiure  20:  Azimuth  Angle  Across  the  Optimized  Path  for  Low  Observability 


Figure  21;  Ground  track  example  5 


Figure  22:  Signature  example  5 
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